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ABSTRACT 

Following an initial explosion that might be launched either by magnetic interactions or neutrinos, 
a rotating magnetar radiating according to the classic dipole formula could power a very luminous 
supernova. While some 56 Ni might be produced in the initial explosion, the peak of the light curve 
in a Type I supernova would not be directly related to its mass. In fact, the peak luminosity would 
be most sensitive to the dipole field strength of the magnetar. The tail of the light curve could 
resemble radioactive decay for some time but, assuming complete trapping of the pulsar emission, 
would eventually be brighter. Depending on the initial explosion energy, both high and moderate 
velocities could accompany a very luminous light curve. 

Subject headings: supernovae: general; hydrodynamics, shock waves, turbulence 



1. INTRODUCTION 



The role of rotation in powering t he explosion of super- 
novae has long been debated (e . g., [A kiyama et al. 120031: 
Hovld 119461 : iLeBlanc fe Wilson! [19701 : lOstriker fc Gunr 



19711 ). Most recent work has focused on the possibility 
that a rotating neutron star could, by way of a mag- 
netic interaction, be the energy source for exploding a 
massive star. While that issue is far from resolved, and 
recent headway has been made in exploding these same 
stars using neutrino transport, it is certain that a large 
fraction of supernova explosions produce rotating neu- 
tron stars and that those neutron stars frequently have 
large magnetic fields. Magnetars are a class of neutron 
stars with field strengths 10 14 to 10 15 Gauss and more. 
They may constitute 10% of the neutron star birthrate 
(|Kouveliotou et all I1998T ) . It is likely that fast rotation 
in the collapsing iro n core is responsible for crea ting the 
large magnetic field ([Duncan fc T hompson 1992), so it is 
reasonable to expect that the birth of rapidly rotating, 
highly magnetic neutron stars is commonplace. 

It is also generally assumed that these rapidly rotating 
neutron stars are magnetically braked by dipole emis- 
sion early on, accounting for the slow periods observed 
in anom alous x-ray pulsars and s o ft gamma-ray re- 
peate rs ([Duncan fc Thompson! 119921 : iKouveliotou et all 
1998). Since the initial rotational energy of such stars at 
birth must have been large and 1000 years later is small, 
where did the difference go? Might it have been emitted 
in some observable form? 

As a rough approximation, assume that the neutron 
star radiates its rotational energy away at a rate given by 
the traditional dipole formula for pulsars. For a typical 
moment of inertia of 10 45 g cm 2 , the rotational energy of 
a neutron star with period, P ms , in milliseconds is 



2 x 10 52 P" 2 erg. 



(1) 



The approximate energy loss for dipole radiation is given 
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by the Larmor formula (e.g.. lLang[ [l980). 



— — = — ^ (BR Sin a) — 

dt 3c 3 v V p 

« 10 49 Bf 5 P" 4 erg s- 1 . 



(2) 



Here B i5 is the surface dipole field in 10 15 Gauss, R » 
10 6 cm is the neutron star radius, and a is the inclination 
angle between the magnetic and rotational axes, taken 
arbitrarily to be 30 degrees. 

At face value, these simple formulae suggest large lu- 
minosities for magnetars during the first few weeks of 
their lives. For example, the magnetar in SGR 1627-41 
is estimated to have a current spin-down luminosity of 
~ 4 x 10 34 erg s~ x at an age of 2.2 k y implying a magneti c 
field BSina - 2 x 10 14 Gauss plsposito et all 120091) . 
Extrapolated back to when the magnetar was 10 days 
old, equation [2] implies a luminosity of 2 x 10 44 erg s _1 . 
Provided the initial rotation rate was rapid compared 
with its value at 10 days, this result is independent of 
the initial rotation rate. 

Here we explore the observational characteristics of 
Type Ib/c supernovae in which pulsars with magnetar- 
like characteristics have been embedded. The bright- 
est supernovae actually come from neutron stars with 
fields that, for a magnetar, are relatively modest, ~ 10 14 
Gauss. Larger fields imply that most of the magnetar's 
rotational energy is dissipated early on, adding to the 
kinetic energy and mixing, but not appreciably affect- 
ing the luminosity. These luminous magnetar-powered 
supernovae might easily be confused with other forms 
of "hypernovae" where the luminosity has a radioactive 
origin. Alternatively, the lack of such emission con- 
strains the existence of any rapidly rotating magnetar. 
This could be an important constraint in the context of 
gamma-ray bursts where the the possibility of a magne- 
tar power source is currently debated. 

2. A MODEL EXPLOSION 

The presupernova model adopted for our first study is 
a 4.37 Mq Wolf-Rayet star derived f rom the evolution 
of sin gle 35 M Q main sequence star ([Wooslev fc Hegerl 
l2007f ). Magnetars may originate fro m a population this 
massive or larger ([Muno et al.l 12006), and the relatively 
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Fig. 1. — Final composition of the explosion whose light curve 
was shown in Fig. [2] The kinetic energy of the explosion is 1 .2 X 10 51 
erg without the pulsar, 1.6 X 10 51 erg with it. 
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small helium core mass will demonstrate the possibility 
of a very luminous supernova from a light progenitor. A 
heavier model is explored later. This star was exploded 
with a piston located at 1.89 Mq, the base of the former 
oxygen burning shell where the entropy per baryon was 
S/NaIs = 4.0. The piston was moved with a speed such 
as to impart a final kinetic energy (without any pulsar 
addition) of 1.2 x 10 51 erg and produce 0.05 Mq of 56 Ni 
(Fig. Q]). The specific nature of this initial explosion is 
left unspecified. It could have been the result of suc- 
cessful neutrino-mediated transport or a consequence of 
early energy deposition related to the formation of the 
magnetar itself. 

Ten seconds after the explosion, when all important 
nucleosynthesis had ceased and the shock wave was still 
about 0.5 Mq beneath the stellar surface, energy depo- 
sition from an assumed embedded magnetar was turned 
on. A magnetic field of 10 14 Gauss and an initial rotation 
period of 4.5 ms were assumed, corresponding to a rota- 
tional kinetic energy of 10 51 erg. Energy from the mag- 
netar was deposited in the inner 10 zones of the model 
at a rate given by equation [2] and complete trapping was 
always assumed. While the form of the energy deposi- 
tion was unspecified, it could be accelerated particles or 
radiation of any wavelength to which the supernova is 
initially opaque. The initial luminosity of the assumed 
pulsar was 10 45 erg s _1 (equation [2|) which declined to 
10 42 erg s _1 after 5.9 x 10^ s. After 2 x 10 11 s the rotation 
period was 1.0 s. 

The resulting light curve is shown in Fig. O At early 
times the effect of the magnetar was negligible. Later, 
some of the energy went into accelerating the ejecta to 
higher speeds and creating a "bubble" inside the super- 
nova. Later the energy was mostly radiated away and, at 
late times, the supernova luminosity tracked equation O 
The light curve peaked at 7.3 x 10 43 erg s -1 (about 5 
typical Type la supernovae) at 37 days. At that time, 
the pulsar was depositing 7.5 x 10 43 erg s _1 . The rise 
time from 10 42 erg s _1 was 24 days. 

Also of some interest is now the magnetar emission 
modifies the dynamics of the explosion. Without the 
pulsar, the explosion energy was 1.2 x 10 51 erg; with it 
the explosion energy is increased to 1.6 x 10 51 erg. That 
is about 40% of the magnetar energy went into kinetic 
energy of expansion, while the remainder went into radi- 
ation. Energy was not deposited uniformly in the ejecta 



Fig. 2. — Light curve of a 4.37 Mq Wolf-Rayct star in which a 
magnetar (10 14 Gauss; E = 10 51 erg) emitting dipole radiation 
is embedded. Complete trapping of the magnetar emission is as- 
sumed at all times. Shown for comparison is the energy released 
by the decay of 4 Mq of 56 Ni. If the light curve at peak had been 
powered by radioactivity it would have required 7 Mq . The second 
solid line, slightly displaced to later times, shows the light curve 
of a heavier star (7.29 Mq) exploded with the same energy and 
containing the same magnetar (see §|2}. 



though. Because the pulsar was a central source of heat- 
ing, its energy went into inflating a bubble of lower den- 
sity material inside the supernova. The boundary of this 
bubble is a dense shell in which resides about a solar 
mass of the ejecta (Fig. [3]). Other calculations, not illus- 
trated here, show that if one assumes a larger magnetar 
rotational energy or stronger magnetic field, the velocity 
and mass of this thin shell are increased. 

In three dimensions, this shell is probably unstable 
and, rather than being ejected as a thin spherical shell, 
the ejecta are probably mixed. This will affect the light 
curve, perhaps making it rise at an earlier time. Multi- 
dimensional studies of the coupled radiation transport 
and hydrodynamics are needed and are feasible, but are 
postponed for now. This mixing will also have important 
implications for the spectrum, especially at late times. 

3. A HEAVIER PROGENITOR AND A STRONGER FIELD 

The same magnetar in supernovae with different 
masses and kinetic energies will give supernovae with dif- 
ferent properties. Consider first the effect of the same 
magnetar as in the previous section (_E ro t = 10 51 erg, 
_£?i5 = 0.1) embedded in a Type lb supernova with mass 
7.29 Mq. The presupernova star here was the remnant 
of a 60 Mq main sequenc e star of solar met a llicity , again 
taken from the study of IWooslev fc Hegeri (|2007l ). The 
star was exploded, as before, with a piston at the base 
of the oxygen shell (1.59 Mq) and given a final kinetic 
energy of 1.2 x 10 51 erg. Explosive nucleosynthesis pro- 
duced 0.11 Mq of 56 Ni though the effect of 56 Ni decay 
was not included in the plots shown. 

Fig. [5] already showed the resulting light curve. It is 
very nearly the same as for the lower mass model. With 
its larger mass and similar kinetic energy, the supernova 
expands slower, but the emerging light is only slightly 
delayed and on the tail, the light curves are identical. 
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Fig. 3. — Final velocity for two explosions in which the pulsar 
energy deposition was included (solid line) or was not (dashed line). 
In one dimension, the energy deposited in the inner ten zones blows 
a bubble resulting in a high density shell moving at nearly constant 
speed. In three dimensions, this shell is probably unstable and 
mixing would occur. 
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Fig. 4. — Light curves for three models based upon the 1.2 xlO 51 
erg explosion of a 7.29 Mq Wolf-Rayet star. The embedded mag- 
netar had a field strength of 5 X 10 14 Gauss and initial rotational 
energies of 2 X 10 50 (solid line) and 5 X 10 50 erg (dotted line) or 
7 X 10 14 Gauss and 5 X 10 50 erg (dot-dashed line). The dashed line 
shows the energy that would result from the decay of 0.3 Mq of 
56 Ni. 



Fig. 2] shows what happens if the magnetar is assumed 
to have a stronger magnetic field, Bi$ = 0.5 and 0.7. 
More of the rotational energy is deposited early on and 
contributes to the expansion rate. Consequently the light 
curve at peak is fainter and it declines more rapidly. 
Smaller initial rotation rates (2 x 10 50 and 5 x 10 50 erg 
corresponding to periods of 10 ms and 6.3 ms) were em- 
ployed, but the answer is not very sensitive to that. The 
extra energy would just accelerate the ejecta. Since the 
resulting thin spherical shells are probably not physical, 
they are not considered here. 

The light curves in Fig. 2] resemble, superficially at 
least, t hat of SN 1998b w which accompanied GRB 
980425 (jPatat et al J 1200 If ). The long tails that in fact 
go like t~ 2 , parallel for some time that expected from 
the decay of 56 Co and if the light curve for the first year 



was all one had to go by, one might easily confuse a 
pulsar-powered model for one powered by radioactivity. 
Since the explosion energy and mechanism is left unspec- 
ified here, one could invoke a much higher kinetic energy 
(accompanying the birth of the magnetar?) and obtain 
spectra like that seen in SN 1998bw. The light curve 
at peak would then be powered by a combination of ra- 
dioactivity and pulsar emission. 

However, th e luminosity of SN 19 98bw at times later 
than one year (jSollerman et all2 002) argues against such 
a hybrid model. By 1000 days after the explosion, the 
supernova had declined to ~ 10 38 erg s _1 , whereas a 
pulsar with field strength 5 x 10 14 Gauss would still be 
radiating 2 x 10 40 erg s _1 . In fact, a strict application of 
equation [2] would imply that the pulsar must have been 
born with a field strength in excess of 8 x 10 15 Gauss in 
order to satisfy the observational limit at late times, but 
then the magnetar would contribute all of its rotational 
energy to the explosion and none to the light curve. 

Of course equation [2] is an approximation whose va- 
lidity can be questioned, especially at early times, and 
the magnetic field strength and orientation need not be 
constant with time. Given the freedom to pick a field 
decay rate, the entire light curve could be fit, but this 
seems somewhat contrived. Fitting the late time light 
curve with a combinati on of 56 Co and 57 Ni de cay seems, 
for now, more natural (jSollerman et a l. 2002). 

4. CONCLUSIONS 

Pulsars with magnetar-like magnetic fields can con- 
tribute appreciably to the light curves of Type lb and 
Ic supernovae. In fact, it would be surprising if they 
never did. For moderate field strengths, the light curve 
is very luminous and long lasting and might be confused 
with those of pair-in stability supe r novae or circumstellar 
interaction. Indeed, iMaeda et al.l (|2007f ) have suggested 
a pulsar as the possible source powering the second (prin- 
cipal) maximum of the light curve of SN 2005bf. 

If magnetars are the central engine that pow- 
ers the long-soft clas s of gamma-ray bursts (e.g., 
IWooslev fc Bloomll2006f ). then it is reasonable to expect 
that they may contribute to the light curves of the su- 
pernovae that accompany them. On the one hand, this 
might facilitate the magnetar paradigm because the pro- 
duction of th e necessary large amoun t of 56 Ni has proven 
problematic (jBucciantini et all [2009). Having an alter- 
nate explanation that involves magnetar energy input 
would solve this problem. On the other hand, if a magne- 
tar contribution to the light curve can be ruled out based 
on the spectrum and late time light curve, one must won- 
der how the magnetar is so effectively concealed. Is the 
rotational energy extracted with such high efficiency in 
the first few seconds that the magnetar forever after- 
wards rotates slowly, or do rotating magnetars not emit 
according to the popular dipole formula during the first 
year? 

Perhaps they do not. The absence of a pulsar contri- 
bution to the luminosity of typical supernovae is easy to 
understand. A pulsar born with a 14 ms period (10 50 
erg), would only have a pulsar luminosity of 6 x 10 erg 
s during the time when it is bright. This is small com- 
pared with the 10 42 - 10 43 erg s _1 resulting from shock 
energy released by recombination (Type II supernova) 
or radioactivity (Type I supernova). But the extremely 
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faint optical luminosity of SN 1987A, < 8 x 10 33 erg s -1 
seventeen years after its birth (|Graves "et~aLl[2005l) is very 
difficult to reconcile with any model with a young active 
pulsar. While dust extinction could be considerable, the 
lack of a point source in SN 1987A remains a mystery. 
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